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Abstract 

In conservation genetics and management, it is important to understand the contri- 
bution of historical and contemporary processes to geographic patterns of genetic 
structure in order to characterize and preserve diversity. As part of a 10-year mon- 
itoring program by the Government of Newfoundland and Labrador, Canada, we 
measured the population genetic structure of the world's most northern native pop- 
ulations of brook trout (Salvelinus fontinalis) in Labrador to gather baseline data 
to facilitate monitoring of future impacts of the recentiy opened Trans-Labrador 
Highway. Six-locus microsatellite profiles were obtained from 1130 fish representing 
32 populations from six local regions. Genetic diversity in brook trout populations 
in Labrador (average He = 0.620) is within the spectrum of variability found in 
other brook trout across their northeastern range, with limited ongoing gene flow 
occurring between populations (average pairwise F sj = 0.139). Evidence for some 
contribution of historical processes shaping genetic structure was inferred from an 
isolation-by-distance analysis, while dual routes of post-Wisconsinan recoloniza- 
tion were indicated by STRUCTURE analysis: K = 2 was the most likely number of 
genetic groups, revealing a separation between northern and west-central Labrador 
from all remaining populations. Our results represent the first data from the nuclear 
genome of brook trout in Labrador and emphasize the usefulness of microsatellite 
data for revealing the extent to which genetic structure is shaped by both historical 
and contemporary processes. 



Introduction 

The extent and structure of genetic variation across a species' 
range are informative about both historical processes, such 
as the pattern of recolonization after previous glacial max- 
ima ( Avise 2000) , and contemporary ones, including ongoing 
dispersal and gene flow among populations. From a conser- 
vation and management perspective, it is important to under- 
stand the contribution of each type of process to the observed 
geographic pattern of genetic structure in order to identify 



evolutionarily distinct lineages and maximally preserve di- 
versity. Effective management requires the maintenance of 
genetic diversity to sustain a population's ability to evolve 
(Frankham et al. 2002) and persevere through environmental 
changes (Moritz and Faith 1998; Avise 2004), such as climate 
change, biological invasions, or anthropogenic threats like 
pollution, habitat fragmentation, or over-exploitation. 

Brook trout (Salvelinus fontinalis Mitchill, 1814; Fig. 1) 
are primarily freshwater salmonid fish (Scott and Crossman 
1 973) that inhabit the clear, cool waters of streams, rivers, and 
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Figure 1. The study species, brook trout (Salvelinus fontinalis). Photo- 
graph by Donald Keefe. 



lakes (Behnke 1972). Coastal populations with ocean access 
are sometimes characterized by anadromous forms of this 
species (Power 1980), contributing to its ecological complex- 
ity. Due to their economic value as game fish, brook trout 
have been widely introduced to temperate regions outside 
their native range and now exist as numerous naturalized 
populations. Endemic to eastern North America (McCrim- 
mon and Campbell 1969), the native range of the species 
includes much of eastern Canada, and extends southward to 
the northern tip of Georgia. Labrador, the mainland com- 
ponent of the province of Newfoundland and Labrador in 
the Atlantic region of Canada, constitutes the northeastern- 
most portion of the native range. Labrador comprises the 
easternmost part of Canadian Shield and consists of two cli- 
mate types, polar tundra in the north and subarctic in the 
south. 

In Labrador, brook trout are the most important game 
fish (DFO 2007), and are under increased anthropogenic 
threat due to the recent opening of the Trans-Labrador High- 
way (TLH). This threat consists of two parts. First, culvert 
placement along the expanse of the highway has been poor 
(Gibson etal. 2005), and the habitat fragmentation that the 
highway imposes could interfere with migration, thus re- 
ducing gene flow and subdividing current populations into 
smaller isolated populations. Previous studies have reported 
the loss of brook trout from much of their native range 
(50% of brook trout have disappeared from their histori- 
cal habitat in Connecticut alone) (Hudy etal. 2008) likely 
due to land development and concomitant habitat alteration 
(Stranko etal. 2008; Kanno etal. 2010; Waco and Taylor 
2010). Second, fishing pressure may increase as lakes be- 
come more easily accessible; previously, the majority of the 
lakes in Labrador were accessible by float plane only. The 
Newfoundland and Labrador Department of Environment 
and Conservation has established a 10-year monitoring pro- 
gram of brook trout populations in proximity of the newly 
developed highway in Labrador. Assessment of the current 
population genetic structure throughout Labrador is crucial 
in order to lay the groundwork for future temporal sam- 
ples of the same populations for monitoring of population 
trends. 
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Several studies have established the role of contemporary 
factors in shaping the genetic structure of brook trout pop- 
ulations across portions of the species' native range, includ- 
ing Maine (Castric etal. 2001), Quebec (Angers etal. 1995, 
1999; Angers and Bernatchez 1998; Fraser etal. 2004), New 
Brunswick (Rogers and Curry 2004), and insular Newfound- 
land (Adams and Hutchings 2003). Ecological variables were 
found to affect the patterning of genetic structure in cer- 
tain cases. For example, in Miramichi River, New Brunswick, 
population structure was dictated by habitat selection and the 
life history strategy of anadromy (Rogers and Curry 2004), 
while the structure of brook trout populations in Mistassini 
Lake, Quebec, was attributed to dispersal patterns, specifically 
male-biased dispersal among tributaries and female-biased 
dispersal from tributaries to outflow populations (Fraser 
etal. 2004). Another major factor that shapes the popula- 
tion genetic structure of freshwater fish such as brook trout 
is the structure of the habitat. In particular, the connectiv- 
ity of lakes within watersheds has been found to influence 
the pattern of genetic diversity of brook trout on a mi- 
crogeographic scale (Angers etal. 1995, 1999; Angers and 
Bernatchez 1998; Castric etal. 2001; Adams and Hutchings 
2003). 

A species' genetic structure is also strongly influenced 
by historical events like the Wisconsinan glaciation which 
reached its maximum 18,000 years ago. Relative to recently 
reestablished northern populations, populations that per- 
sisted in isolation in southern refugia have had time to ac- 
cumulate variation (Hewitt 1996, 1999; Bennett 1997) while 
simultaneously diverging from populations in other refu- 
gia (Hewitt 1996, 2004). This signature is reflected in the 
distribution and genetic structure of populations today - 
northern regions typically have fewer species (Pielou 1991) 
characterized by lower genetic diversity, a pattern that is re- 
ferred to by Hewitt (1996) as "southern richness, northern 
purity." The pattern of population genetic structure there- 
fore reveals the direction of colonization and the contribu- 
tion of distinct refugial origins to the present-day species' 
distribution. 

The pattern of colonization of brook trout throughout 
their native range was assessed by Danzmann et al. (1998) on 
the basis of mtDNA diversity. Brook trout in eastern Canada 
were shown to possess very little genetic variation relative 
to the rest of the range, comprising only 13 haplotypes of a 
possible 61. This part of the range was primarily colonized 
by haplotype " 1 ," postulated to have been present in both the 
Mississippian and Atlantic refugia. However, the presence 
of 10 private haplotypes combined with the lack of haplo- 
type "2" fish suggests that fish in this region may instead 
have originated from a more northeasterly refugium (Danz- 
mann etal. 1998), also referred to as the Acadian refugium 
(Schmidt 1986). Although Danzmann etal. (1998) did not 
consider Labrador brook trout, Black etal. (1986) inferred 
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from species distributional patterns that brook trout recol- 
onized inland Labrador via Quebec from a combination of 
the Mississippian and Atlantic refugia while simultaneously 
recolonizing coastal regions of Labrador from the Atlantic 
refugium. This pattern is inconsistent with Danzmann et al. 
(1998). 

As part of the Department of Environment and Conser- 
vation's monitoring program, we sought to determine the 
genetic structure of brook trout populations in Labrador as 
recorded by nuclear microsatellite loci. Such analyses will 
provide the baseline information about current population 
structure needed to monitor population trends impacted by 
the TLH. In addition, we test whether patterns of genetic 
variation are influenced more strongly by historical or con- 
temporary events. Since Labrador brook trout populations 
are at their northern range extreme, we hypothesize that they 
maintain the genetic signature of postglacial recolonization. 
Since the mtDNA phylogeography and distributional pat- 
terns of fishes in Labrador did not clearly elucidate postglacial 
colonization routes of brook trout, we also investigated the 
number and source of refugial populations that colonized 
Labrador. 

Methods and Materials 
Sample collection 

Brook trout caudal fin clippings from fish in 31 lakes within 
18 watersheds from six regions across eastern Canada (north- 
ern, west-central and southeastern Labrador, insular New- 
foundland, Nova Scotia, and New Brunswick; Fig. 2) were 
sampled via gill netting by provincial wildlife biologists from 
2003 to 2008 (1094 samples). Fin clippings were stored at 
-20° C. DNA extractions of fish from Atikonak Lake in west- 
central Labrador (36) were provided by Dr. Steven M. Carr 
(Memorial University of Newfoundland and Labrador). 

Microsatellite analysis 

Genomic DNA was extracted from ~25 mm 2 of the cau- 
dal fin clipping using the Qiagen QIAmp DNA Mini Kit 
(Qiagen Inc., Mississauga, Canada) following the tissue pro- 
tocol. Nine microsatellite loci were initially tested: S/ol2, 
S/o 18 (Angers etal. 1995), MST85 (Presa and Guyomard 
1996), S/oD91, S/oB52, S/oC129, S/oDlOO, S/0C86, S/oD75 
(T.L. King, US Geological Survey, unpubl. data). PCRs were 
carried out in an Applied Biosystems GeneAmp PCR System 
9700 Thermal Cycler (Applied Biosystems Inc., Carlsbad, CA) 
and contained IX PCR Master Mix (Promega Corp., Madi- 
son, WI), 2 /xL 10 mM of each primer and 2 /xL genomic 
DNA (2-570 ng). Four fluorescent dyes were used to label 
forward primers as follows: PET - Sfoll, S/oDlOO; VIC - 
S/0I8, MST85, and S/oC129; 6-FAM - S/oD91, S/oB52, and 
S/0C86; NED - S/oD75. For primer pairs Sfol2, SfolS, and 
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MST85 the PCR profile was as follows: 95°C for 4 min; 35 
cycles of 95°C for 1 min, 55°C for 45 s, and 72°C for 45 s; 
72° C for 7 min. For the remaining loci the PCR profile was: 
95°C for 4 min; 35 cycles of 94°C for 45 s, 56°C for 45 s and 
72°C for 1.5 min; 72°C for 5 min. Prior to electrophoresis, 
reactions were diluted (1:4), and 1 /xL was added to 8.8 /xL 
of formamide and 0.2 /xL of LIZ-500 size standard. Loci were 
multiplexed as follows: (1) S/o 12, S/0I8, S/oD91; (2) S/oB52, 
S/oC129, S/oDlOO; (3) S/0C86, S/oD75,MST85. Samples were 
electrophoresed in an Applied Biosystems 3730 DNA An- 
alyzer using GeneScan Software (Applied Biosystems Inc., 
Carlsbad, CA). The six loci with the greatest amplification 
success rate (S/0I8, S/oB52, S/oC129, S/oDlOO, S/0C86, and 
S/oD75) were selected for subsequent analysis. 

Data analysis 

Diversity at each locus was quantified by calculating observed 
(H 0 ) and expected (H E ) heterozygosity and number of alle- 
les (N A ) in FSTAT 2.9.3 (Goudet 1995). The inbreeding co- 
efficient (Fis) at each locus was estimated using the Weir and 
Cockerham (1984) algorithm as implemented in FSTAT (120 
randomizations). Linkage disequilibrium among all pairs of 
loci was evaluated with the log likelihood ratio test (G-test) 
of GENEPOP 1.2 (Raymond and Rousset 1995) using the 
default Markov chain parameters. Diversity measures (N A , 
Ho, H E ) and allelic richness (A) were calculated in FSTAT 
for each population (defined as all fish sampled in a lake), as 
well as each region (defined by Anderson [1985] according to 
geology, vegetation, and fish distributions, with watersheds 
in each region draining to a specific inlet or region of coast). 
The rank-based Kruskal-Wallis test was used to determine if 
differences between overall H E and A of each lake were sig- 
nificant. FIS values for each locus in each lake were calculated 
in FSTAT. 

To determine whether mutation as well as drift has influ- 
enced differentiation, we used the allele size randomization 
test (following Hardy etal. 2003) ofSPAGeDi 1.3 (Hardy and 
Vekemans 2002). Randomly permuting alleles (2000 itera- 
tions) provides a 95% confidence interval of the simulated 
distribution of Rst values (or pRST). If _R S t > fsT> Rst: will 
fall outside of the pRsj 95% confidence interval, and differ- 
entiation has been influenced by mutation. If it falls within 
this range, _R S t does not significantly differ from F S j, indicat- 
ing that mainly drift contributes to differentiation. Since we 
found that _R ST (0.153) was not significantly different from 
Fst (-Fst = 0.144; P = 0.115) (Table 3), pairwise _F S t values, 
as estimated in ARLEQUIN 3.1 (Excoffier etal. 2005), were 
used in subsequent analyses. Due to greater variance Rst is 
less reliable at detecting population differentiation than F ST 
(Sefcetal. 2007). 

To assess the partitioning of genetic diversity among pop- 
ulations, an analysis of molecular variance (AMOVA) was 
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Figure 2. Lakes sampled for brook trout across eastern Canada. Population codes are given in Table 2. 



performed in ARLEQUIN 3.1 (Excoffier etal. 2005), with 
three hierarchical levels to quantify the variance among 
(1) lakes within regions; (2) watersheds within regions; 
and (3) lakes within watersheds. In addition, Nei etal.'s 
(1983) D A distance (calculated with MICROSATELLITE 
ANALYSER 4.05 [Dieringer and Schlotterer 2003]) was used 
to assess distance-based relationships among watersheds. 
The neighbor-joining dendrogram was generated using the 
NEIGHBOUR and CONSENSE modules in PHYLIP 3.69 
(Felsenstein 1989), and visualized in TREEVIEW 1.4 (Page 
1996). 

To characterize the correlation of genetic distances and ge- 
ographical distances, we conducted an isolation-by-distance 
analysis. Pairwise estimates of Fst/(1-Fst) were regressed 
against geographical distance (Rousset 1997) using the 
ISOLDE program within GENEPOP. Distances were mea- 
sured in two ways: (1) in a straight line overland between 
lakes, and (2) between lakes along extant waterways, to deter- 
mine if historical or contemporary influences play the more 
significant role in shaping observed genetic structure. Dis- 
tances were measured using ESRI Arc View 9 (Redlands, CA). 

STRUCTURE 2.1 (Pritchard etal. 2000) was used to iden- 
tify the number of genetically distinct clusters (K) across 



all sampled lakes. First we wanted to discern between the 
possible routes of postglacial recolonization, either from a 
combination of Atlantic and Mississippian refugia, or from 
the Acadian refugium. Secondly we wanted to evaluate the 
overlay of contemporary patterns of gene flow. For each value 
of K, five iterations were run to assess convergence of the like- 
lihood with a burn-in period of 50,000, followed by 200,000 
iterations for values of K = 1 through 20. Each simulation was 
performed using an ancestry model incorporating admixture 
and correlated allele frequencies, without prior population 
information (Falush et al. 2003). K was selected in two ways: 
(1) as the highest value of In Pr(XI^T), as recommended by 
Pritchard et al. (2000) and (2) according to an ad hoc statistic 
described by Evanno et al. (2005). 

Results 

Per-locus diversity and linkage 
disequilibrium 

A total of 1130 brook trout from 32 lakes within 19 water- 
sheds across six regions of eastern Canada were genotyped at 
six microsatellite loci. An average of 17 alleles (N A ) across all 
six loci was found, with a range of eight alleles at S/oC129 to 
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Table 1. Allelic variation at each microsatellite locus. Number of alleles 
(W A ), observed (H 0 ) and expected (H E ) heterozygosity and F B (Wright's 
inbreeding coefficient). Significant heterozygote excess/deficit at P < 
0.05 after Bonferroni correction (a = 0.00833) denoted with an asterisk. 



Locus 


H 0 


H E 


W A 


Fis 


SMS 


0.757 


0.597 


19 


0.211* 


SfoB52 


0.643 


0.576 


15 


0.105* 


SfoC129 


0.527 


0.528 


8 


-0.001 


5foD75 


0.698 


0.616 


17 


0.118* 


5foD91 


0.872 


0.732 


31 


0.161* 


SfoD100 


0.779 


0.807 


12 


-0.036 


Weighted averages 


0.717 


0.642 


17 


0.098* 



31 alleles at S/oD91 (Table 1). Observed heterozygosity (H 0 ) 
ranged from 0.528 (S/oC129) to 0.807 (S/oDlOO), with an av- 
erage of 0.642. A preliminary indication of genetic structure 
was detected by deviations from Hardy- Weinberg equilib- 
rium. Fi S was significantly greater than zero for four loci 
(S/ol8, S/oB52, S/bD75, SfoD91) and overall after correction 
for multiple comparisons (P = 0) (Table 1). Significant link- 
age disequilibrium after Bonferroni correction was found for 
only three of 480 pairwise comparisons of loci by the log 
likelihood ratio test (G-test) of GENEPOP (P < 0.0001). 

Within population variation 

Levels of genetic variation were highest in Farnham Brook in 
Nova Scotia (H E = 0.727, A = 4.910), and lowest in Dead 
Dog Pond in central Labrador (H E = 0.370, A = 2.616) 
(Table 2). When populations were assigned to their corre- 
sponding region, genetic diversity was highest in Newfound- 
land in terms of expected heterozygosity (H E = 0.702) and 
in New Brunswick based on allelic richness (A = 8.347) 
and lowest in northern Labrador (H E = 0.682, A = 6.272) 
(Table 2). The Kruskal-Wallis analysis of variance showed 
no significant difference between average expected heterozy- 
gosity (over all loci) and allelic richness on either a per-lake 
or a regional basis (P > 0.05). However, there was a signif- 
icant negative correlation between allelic richness and lati- 
tude (R 2 = 0.214; P = 0.009), as well as a significant positive 
correlation between heterozygosity and longitude in brook 
trout populations found in west-central and southeastern 
Labrador (R 2 = 0.2476; P = 0.0184), after outliers were 
removed. No significant deviations from Hardy-Weinberg 
equilibrium were observed per locus among lakes or overall 
(P> 0.05). 

Population genetic structure 

The mean F ST measured among pairs of lakes was 0.139 and 
comparisons were significantly (P < 0.05) positive except 
those involving Dead Dog Pond (which had the smallest sam- 
ple size of N = 11). The significant pairwise Fst measures 



ranged from 0.014 between Gilbert Lake and Feeder Pond 
(both in southeastern Labrador), to 0.374 between Atikonak 
Lake (west-central Labrador) and Alder Brook (Nova Sco- 
tia). The average pairwise F S j among lakes in a watershed 
measured 0.070. The isolation-by-distance analysis among 
all populations revealed a significant correlation between 
Fst/Q-Fst) and both overland and current waterway dis- 
tances, but was more highly significant for overland distances 
(P = 0 vs. P = 0.002). 

Both AMOVA analyses detected significant geographic 
structuring within and among populations and regions (P < 
0.05; Table 4). The first analysis defined groups by lakes, then 
watersheds. Most of the variation was detected within lakes, 
followed by among watersheds, then among lakes within wa- 
tersheds. In the second analysis, fish were grouped by water- 
shed, then into regions. Most of the variation was detected 
within watersheds, followed by among watersheds within re- 
gions, and the least was among regions. Overall, the majority 
of variation is within lakes, and the next most significant fac- 
tor explaining variation is watershed. There were very low 
but significant levels of differentiation (P = 0.03) among 
regions. 

Genetic structure was also investigated using Nei's genetic 
distance (D A ) (Nei etal. 1983), the most appropriate dis- 
tance measure when drift is the main force contributing to 
differentiation (Paetkau etal. 1997). Following the AMOVA 
results, fish were grouped into watersheds for this analy- 
sis. The neighbor-joining dendrogram of D A distances illus- 
trates overall low genetic differentiation among watersheds 
with some regional pattern (Fig. 3). All watersheds in north- 
ern Labrador constitute one cluster. Watersheds from New 
Brunswick, Nova Scotia, and Newfoundland, with the ex- 
ception of Middle Brook, cluster together. All watersheds 
in southeastern Labrador that drain directly into the At- 
lantic Ocean (south of Sandwich Bay) cluster together. Water- 
sheds in southeastern and west-central Labrador that drain 
from Sandwich Bay and Lake Melville, respectively, cluster 
together. 

From the STRUCTURE analysis, under the Pritchard 
etal. (2000) criterion, the maximum value of In Pr(XI^C) 
(-19197.2) corresponds to K = 14 (Fig. 4). Upon mapping 
the bar plots of individual estimated cluster membership co- 
efficients of K = 14 (Fig. 5A), lakes within the same wa- 
tershed tend to have similar cluster assignments throughout 
west-central and southeastern Labrador. Exceptions to this 
pattern include Nikki's Pond (Kenamu River), which has a 
cluster assignment distinct from any other lake or watershed; 
The Right Lake (Traverspine River), which clusters with Ea- 
gle River; Nippard's Lake (Eagle River); and No Boat Pond 
(Traverspine River), which have similar cluster assignments 
distinct from their source watersheds. We note that one clus- 
ter assignment dominates all three watersheds in northern 
Labrador, and the two watersheds in Newfoundland have 
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Table 2. Locations with population code, sample sizes (N), observed (H 0 ) and expected (H E ) heterozygosity, and allelic richness (A, corrected to n = 7) 
of each lake and across each region. 



Region 


Watershed 


Lake 


N 


We 


H 0 


A 


Northern Labrador 




36 


0.654 


0.682 


6.272 


1 


Saputit Brook 


Saputit Lake 


7 


0.714 


0.658 


4.000 


2 


Kogaluk River 


Cabot Lake 


16 


0.683 


0.609 


3.329 


3 


I lad I ivi k Brook 


Walkabout Lake 


13 


0.585 


0.639 


3.715 


West-central Labrador 




207 


0.585 


0.700 


6.985 


4 


Atikonak Lake 


Atikonak River 


22 


0.574 


0.649 


4.245 


5 


Traverspine River 


No Boat Pond 


45 


0.385 


0.470 


3.378 


6 


Traverspine River 


The Right Lake 


17 


0.627 


0.608 


4.013 


7 


Kenamu River 


Mercier Lake 


15 


0.491 


0.616 


3.739 


8 


Kenamu River 


Brennan Lake 


58 


0.722 


0.694 


4.265 


9 


Kenamu River 


Nikki's Pond 


50 


0.623 


0.588 


3.635 


Southeastern Labrador 




734 


0.650 


0.692 


8.1 18 


10 


Eagle River 


Fred's Lake 


53 


0.727 


0.691 


4.491 


11 


Eagle River 


No Name Lake 


27 


0.595 


0.616 


3.935 


12 


Eagle River 


Nap Pond 


33 


0.613 


0.630 


4.238 


13 


Eagle River 


Nippard's Lake 


60 


0.584 


0.568 


3.348 


14 


Eagle River 


Osprey Lake 


37 


0.653 


0.610 


3.542 


15 


Eagle River 


Dead Dog Pond 


1 1 


0.268 


0.370 


2.616 


16 


St. Augustine River 


St. Augustine 


54 


0.586 


0.603 


3.776 


17 


St. Augustine River 


Bog Lake 


59 


0.558 


0.574 


3.633 


18 


Paradise River 


Keith's Lake 


41 


0.603 


0.591 


4.072 


19 


Paradise River 


Crooked Lake 


55 


0.521 


0.591 


4.177 


20 


Alexis River 


Alexis Pond 


45 


0.806 


0.667 


4.056 


21 


Alexis River 


Handkerchief Pond 


29 


0.673 


0.575 


3.839 


22 


Alexis River 


Feeder Pond 


24 


0.627 


0.603 


4.004 


23 


Gilbert River 


Tilt Pond 


46 


0.706 


0.635 


4.450 


24 


Gilbert River 


Gilbert Lake 


60 


0.672 


0.592 


4.048 


25 


St. Lewis River 


Curl's Pond 


61 


0.817 


0.688 


4.457 


26 


St. Mary's River 


Mary's Harbour Big Pond 


39 


0.719 


0.680 


4.632 


Newfoundland 






54 


0.654 


0.702 


6.671 


27 


Middle Brook 


Butt's Pond 


25 


0.639 


0.637 


4.330 


28 


Salmonier River 


Little Gull Pond 


29 


0.662 


0.628 


3.590 


Nova Scotia 






50 


0.690 


0.694 


7.889 


29 


River Denys 


Alder Brook 


30 


0.643 


0.620 


4.391 


30 


Salmon River 


Farnham Brook 


20 


0.769 


0.727 


4.910 


New Brunswick 






49 


0.706 


0.690 


8.347 


31 


Miramichi River 


Moose Lake 


25 


0.647 


0.649 


4.641 


32 


Inner Bay of Fundy 


Walton Lake 


24 


0.767 


0.663 


4.159 



distinct cluster assignments from each other and all other 
watersheds, while Nova Scotia and New Brunswick are char- 
acterized by a heterogeneous cluster assignment distinct from 
the other regions. 

Under the ad hoc criterion of Evanno etal. (2005) K = 2. 
Evanno etal. (2005) recommended this statistic when In 
Fr(X\K) values plateau once the optimal K has been reached, 
as was the case here (Fig. 4). Bar plots of estimated individual 
cluster membership coefficients for K = 2 are mapped in 
Figure 5B. The two clusters consist of (1) Northern Labrador 
and all populations (except Atikonak Lake) that drain into 
Lake Melville and all populations that drain into Sandwich 



Table 3. Summary of allele size permutation test as implemented in 
SPAGeDi showing estimates of R ST , simulated R ST (p/? ST ) values and their 
95% confidence interval and f S T following 2000 allele permutations. 



Locus 


/?ST 


PRST 


95% CI 


F ST 


P 


Sfo18 


0.166 


0.127 


0.045-0.267 


0.156 


0.297 


SfoB52 


0.131 


0.114 


0.054-0.204 


0.130 


0.193 


5foC129 


0.145 


0.150 


0.084-0.200 


0.167 


0.709 


SfbD75 


0.137 


0.152 


0.062-0.267 


0.185 


0.837 


SfoD91 


0.112 


0.076 


0.045-0.166 


0.089 


0.100 


SfoD100 


0.128 


0.088 


0.044-0.169 


0.101 


0.169 


All 


0.153 


0.130 


0.087-0.176 


0.144 


0.115 
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Figure 3. Unrooted neighbor-joining 
denrogram using Nei et al.'s (D A ) distance 
(1 983) relating the 33 brook trout populations 
sampled across eastern Canada. 



Gilbert 
River 




Middle 
Brook 



Eagle 
River 



St. Augustine 
River 



Traverspine 
River 



Atikonak 
River 



Table 4. Hierarchical analysis of molecular variance in brook trout populations across northeastern Canada performed by grouping fish into populations 
(either by lake or watershed), then subsequently into either watersheds or regions. Percentage of total variance (%) and F-statistics (F C t, F sc , Fst) for 
each hierarchical level represented. 





Among 

regions/watersheds 


Among populations 
within 

watersheds/regions 


Within populations 




Comparison 


% F a 


% Fsc 


% 


F ST 


i. Lakes grouped into watersheds 

ii. Watersheds grouped regionally 


7.516 0.075 
1.313 0.013 


6.550 0.071 
10.455 0.106 


85.935 
88.232 


0.141 
0.118 



Bay and (2) the coastal populations of southeastern Labrador, 
insular Newfoundland, Nova Scotia and New Brunswick, and 
Atikonak Lake. 

Discussion 

We investigated population genetic structure of brook trout 
in Labrador and the broader region of northeastern North 
America, including insular Newfoundland, Nova Scotia, and 
New Brunswick, by measuring variation at six microsatellite 
loci. We have collected the first genetic data from the nu- 
clear genome on brook trout in Labrador, providing us with 
baseline data about current population structure needed to 
monitor future impacts of the recently opened TLH. We also 
make inferences about patterns of post-Wisconsinan colo- 



nization of Labrador by brook trout, which are significant 
in the broader context of postglacial recolonization and the 
biogeography of fishes in Atlantic Canada. 

Patterns of intrapopulation genetic 
variability reflect putative postglacial 
recolonization routes 

Levels of expected heterozygosity in brook trout populations 
in Labrador, Canada (average H E = 0.620), are within the 
spectrum of variability observed throughout the northeast- 
ern range (H E = 0.380-0.778 [Angers and Bernatchez 1998; 
Adams and Hutchings 2003; Rogers and Curry 2004; Pois- 
sant etal. 2005]). Since Labrador is located at the northern- 
most edge of the brook trout's native range, diversity was 
expected to be reduced due to relatively recent bottlenecks 
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Figure 4. A plot of the log likelihood of the data (In Pr(XI<» values obtained from the STRUCTURE analysis for /C = 1—15. 



and founder events. However the range of H E reported in 
other studies is so broad and considers so many different 
geographic scales that it is difficult to directly accept or re- 
ject this hypothesis. However, the comparatively moderate 
levels of microsatellite genetic diversity we report is inconsis- 
tent with the only other study to investigate genetic variation 
in brook trout populations in Labrador. In an assessment 
of the intraspecific population genetic structure of Atikonak 
Lake in western Labrador (S. Carr, unpubl. report) only one 
mtDNA haplotype (401 bp of the cytochrome b gene) was 
found among 23 individuals. Similarly a wide-scale survey of 
brook trout RFLP variation showed that 98% of the fish that 
colonized northeastern North America were of a single hap- 
lotype (Danzmann etal. 1998). Clearly, the rapidly evolving 
nuclear microsatellite loci are potentially more informative 
than mtDNA for investigating population structure in this 
recently colonized region. 

As measured by allelic richness, brook trout in northeast- 
ern Canada demonstrate the characteristic "southern rich- 
ness, northern purity" (Hewitt 1996) pattern; A is highest 



in Nova Scotia and lowest in northern Labrador (Table 1). 
In fact there is a significantly negative correlation between 
allelic richness and latitude (P = 0.018); this supports a route 
of postglacial recolonization through northeastern Canada 
from more southern refugia, resulting in fewer alleles in 
northern populations due to a bottleneck effect. 

The negative correlation between expected heterozygosity 
and latitude was not significant (P > 0.05), probably be- 
cause H E is lowest overall in west-central rather than north- 
ern Labrador. Part of west-central Labrador is elevated on a 
plateau, a type of geography with a number of consequences. 
Areas of higher elevations remain glaciated for longer pe- 
riods, and later population settlement means less time to 
accumulate variation. Higher elevations may also have been 
more challenging to colonize, resulting in a lower number 
of founders and smaller effective population size. In ad- 
dition, populations in higher altitudes are more likely to 
be isolated by waterfalls, with net emigration downstream 
again reducing population size and increasing isolation. Re- 
duced heterozygosity at higher altitudes has been reported 
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I 



Nova Scolia 




Figure 5. Geographical distribution of genetic clusters inferred from the STRUCTURE analysis based on K = 14 (A) and K = 2 (B) for each watershed 
sampled. Each bar represents an individual fish with each color representing a genetic group or cluster. Bar charts are sized according to sample size, 
and lakes within the same watershed appear together. 



in a number of fish species (Trinidadian guppies [Poecilia 
reticulata; Shaw etal. 1991]; eastern mosquitofish, [Gambu- 
sia holbrooki; Hernandez-Martich and Smith 1990]; brown 
trout [Salmo trutta; Hamilton etal. 1989]) including brook 
trout in Quebec (Angers etal. 1999). 

Notably, a significant positive correlation was found be- 
tween expected heterozygosity and longitude in brook trout 
populations in west-central and southeastern Labrador. 
Brook trout populations along the eastern coast of Labrador 
whose watersheds drain directly to the Atlantic Ocean (Alexis 
River, Gilbert River, St. Lewis River, St. Mary's River) tend to 
have higher heterozygosities. Previous studies have shown 
that anadromous fish have higher levels of heterozygosity 
than freshwater resident fish (Gyllenstein 1985; Ward etal. 
1994; DeWoody and Avise 2000; Tonteri et al. 2007), a trend 
which is most often attributed to larger long-term effective 
population sizes. Therefore it is possible that brook trout 
populations found in watersheds in close proximity to the 
Atlantic coast are anadromous; a reasonable conjecture given 
that brook trout are known to have alternate life history styles 
when found in coastal areas (Curry etal. 2010). In support 



of this, if brook trout in southeastern Labrador did recolo- 
nize the region from an Acadian refugium (as suggested by 
Danzmann et al. 1998) fish would have had to be euryhaline 
(a trait associated with anadromy) to disperse back to the 
mainland. Southeast coast populations would also have been 
the first recolonized from an Acadian refugium, and hence 
should be the most variable. 

Contemporary patterns of genetic variation 
suggest limited ongoing gene flow 

Pairwise F$j measures in Labrador revealed similar levels 
of differentiation (mean F S j = 0.139 among lakes; data 
not shown) to other studies of brook trout populations, 
including those in Maine, Quebec, and Atlantic Canada 
(f ST = 0.099-0.373; Angers and Bernatchez 1998; Hebert 
etal. 2000; Castric etal. 2001; Castric and Bernatchez 2003). 
The high level of divergence among populations in Labrador 
suggests that there are limitations to ongoing gene flow 
and the hierarchical AMOVA indicates that watershed is 
the predominant level of differentiation (lakes grouped into 
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watersheds: F C t = 0.075; watersheds grouped into regions: 
_F SC = 0.106). Consistent with this, the average pairwise Pst 
among lakes within watersheds was much lower (P S t = 0.070) 
than among lakes overall. This pattern can be explained by 
the freshwater resident lifestyle of the majority of brook trout 
in Labrador, and is consistent with tagging data indicating 
that brook trout are moving among lakes (R. Perry, pers. 
comm.). By way of contrast, a microsatellite and tagging 
study of brook trout in Indian Bay, insular Newfoundland, 
revealed very little movement of brook trout among lakes 
(Adams and Hutchings 2003 ) . The pattern we observed is typ- 
ical of freshwater fish (Gyllenstein 1985; Currens etal. 1990; 
Carvalho etal. 1991; Ward etal. 1994; Estoupetal. 1998; Huey 
etal. 2010) including brook trout (Angers and Bernatchez 
1998) where genetic structure is highly influenced by drainage 
structure. Brook trout may be more mobile in Labrador than 
insular Newfoundland due to its colder climate and harsher 
winter — migratory behavior would allow brook trout to max- 
imally exploit resources in this type of environment (Power 
2002). 

Support for a strong pattern of differentiation among wa- 
tersheds was provided by the STRUCTURE analysis in which 
assignment of individuals into K = 14 clusters distinguishes 



the majority of watersheds from one another. One exception 
to this pattern is that Nikki's Pond (Kenamu River) has a 
unique cluster assignment, suggesting genetic isolation and 
differentiation of this lake. This could be due to a natural 
barrier (e.g., a waterfall), or to the TLH which crosses the Ke- 
namu River watershed along a river that isolates Nikki's Pond 
from Mercier Lake and Brennan Lake (which are genetically 
similar to one another). A second exception to the pattern is 
that the Traverspine River and Eagle River watersheds share a 
cluster assignment, although each contains a single lake with 
a unique cluster assignment. The genetic similarity of these 
two watersheds may be a historical relict reflecting tempo- 
ral instability of watershed structure; as certain lakes within 
these separate watersheds are only ~100 km apart they may 
have been connected at the time of glacial retreat. Ongoing 
gene flow of anadromous fish among watersheds is an un- 
likely explanation, as the river mouths where the watersheds 
drain are far apart (~300 km). 

We also observed that brook trout populations along the 
eastern coast of Labrador whose watersheds drain directly 
into the Atlantic Ocean have STRUCTURE cluster assign- 
ments that are similar to one another, but distinct from the 
rest of Labrador. These are the same watersheds (Alexis River, 
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Gilbert River, St. Lewis River, St. Mary's River) that tend to 
have higher levels of heterozygosity and are hypothesized to 
contain anadromous brook trout. The connectedness of these 
coastal populations, also supported by low levels of differen- 
tiation (pairwise Fst = 0.066) and hence membership in the 
same neighbor-joining dendrogram cluster (Fig. 3), provides 
further evidence of an anadromous life history. 

Implications for management of brook trout 
in southeastern Labrador 

Previous studies have reported extensive loss of brook trout 
from their native habitat due to land development and habitat 
alteration (Hudy et al. 2008; Stranko et al. 2008; Kanno et al. 
2010; Waco and Taylor 2010). The populations surveyed in 
this study are under the same anthropogenic threats due to 
the proximity of the TLH which crossed 20 watersheds in its 
Phase II section (Red Bay to Cartwright) alone (Jacques Whit- 
ford Environment Limited 1998). Therefore it is important 
to collect baseline data on the genetic diversity and structure 
of populations to make management recommendations and 
monitor populations in the future. Previous management de- 
cisions made by the Wildlife Division of the Newfoundland 
and Labrador Department of Environment and conservation 
relied on Adams and Hutchings (2003) which found that 
brook trout do not migrate among lakes. By way of con- 
trast, microsatellite data combined with preliminary tagging 
data (Robert Perry, pers. comm.) suggest that brook trout 
in Labrador do move among lakes within watersheds, and 
that watersheds represent the unit of reproductive isolation. 
Hence it may be more appropriate to direct management 
decisions at the watershed level in Labrador. 

Gibson etal. (2005) reported that 53% of the culverts as- 
sessed along Phase II of the TLH are not suitable for fish 
passage due to either poor design or installation. This causes 
concern with respect to loss of spawning and rearing habi- 
tats, as well the introduction of barriers to migration. Barriers 
to migration are especially worrisome in light of the micro- 
satellite evidence presented here that brook trout in Labrador 
are in fact moving among lakes. Fortunately, the Gibson et al. 
(2005) study led to the improvement of a number of culverts 
along Phase II, as well as the revision and development of 
better stream crossings for Phase III. Regardless, it is impor- 
tant to continue to monitor populations, as easier access to 
remote watersheds can also lead to fisheries' exploitation. 

Evidence for dual routes of post-Wisconsinan 
colonization 

The pattern of microsatellite diversity in Labrador brook 
trout is clearly influenced by contemporary factors, in par- 
ticular watershed structure. However, a role for historical 
events is also evident, as despite structure at the level of wa- 
tershed, certain watersheds are genetically similar, suggest- 
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ing they were colonized by the same founding population, 
and there is little genetic differentiation among regions. The 
isolation-by-distance analysis also indicated that historical 
waterway distances have been more influential than contem- 
porary ones in shaping the present-day genetic structure of 
brook trout populations in northeastern Canada. Consistent 
with this, brook trout populations in the northernmost por- 
tion of their range are evolutionarily young and likely not yet 
in equilibrium between migration and genetic drift. Current 
patterns of microsatellite variability therefore reflect colo- 
nization history overlaid with developing divergence due to 
ongoing limitations to gene flow. 

Microsatellites have been successfully utilized to infer col- 
onization patterns (Angers and Bernatchez 1998; Taylor and 
McPhail 2000; Lu etal. 2001; Koskinen etal. 2002), and can 
reveal patterns that are undetectable with mtDNA or other 
phylogeographic markers (Flanders etal. 2009). This is the 
case for Labrador brook trout; only one mtDNA haplotype 
was found during a survey of genetic variation of freshwater 
fish in Labrador (S. Carr, unpubl. report) and Danzmann 
etal. (1998) reported the predominance of a single mtDNA 
haplotype in the entire northern portion of the range. 

Here we detected evidence for K = 2 genetic groups from 
the STRUCTURE analysis in which we applied the Evanno 
etal. (2005) ad hoc algorithm; one that includes lakes in 
northern and west-central Labrador, and another encom- 
passing lakes in all other regions. This deep level of genetic dif- 
ferentiation suggests that northern and west-central Labrador 
populations of brook trout may have been colonized from a 
different refugial source than the rest of the northeastern 
range. 

Black et al. (1986) proposed that brook trout from the At- 
lantic and Mississippi refugia dispersed through Quebec then 
into Labrador, while other fish from the Atlantic refugium in- 
vaded coastal regions. Alternatively, Danzmann etal. (1998) 
suggested that brook trout in the northeastern parts of their 
native range may have originated from a more northeast- 
erly Acadian refugium. Our results provide support for a 
combination of the two hypotheses; northern and west- 
central Labrador may have been founded by some com- 
bination of fish from the Atlantic and/or Mississippian 
refugia, while southeastern populations were founded by 
those from an Acadian refugium. In support of this, the 
neighbor-joining dendrogram clusters northern Labrador 
populations with those in west-central Labrador (except 
Atikonak Lake). A third hypothesis is that northern and 
west-central Labrador lakes were colonized from a north- 
ern coastal refugium/nunatak located in the Torngat Moun- 
tains in northern Labrador (Pielou 1991). In this scenario 
central Labrador lakes would have been colonized subse- 
quent to northern ones, accounting for their reduced diver- 
sity. Evidence for a Labrador refuge is also supported by the 
phylogeographic origins of Arctic charr (Wilson et al. 1996), 
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as fish in this region were not represented by the Lauren- 
tian lineage which colonized Quebec, New Brunswick, and 
New England. 

Conclusions 

We note that microsatellite variation is sufficient to detect the 
presence of two refugial sources of brook trout in Labrador, 
whereas previous mtDNA studies revealed little to no varia- 
tion (S. Carr, pers. comm.; Danzmann etal. 1998). In addi- 
tion to the strong role of the historical process of recoloniza- 
tion on shaping the pattern of present-day genetic structure 
of brook trout populations in Labrador, we find evidence 
that a contemporary factor, watershed structure, is also im- 
portant. Wildlife management practices relying on the data 
from Newfoundland brook trout have assumed that there is 
no movement of fish among lakes in Labrador; we found 
evidence in contradiction of this which has important im- 
plications for future management decisions with respect to 
the recent opening of the TLH. We point to the importance 
of documenting the genetic structure of a species across its 
entire range rather than extending conclusions from a por- 
tion of the range, especially in species like brook trout that 
are ecologically complex and display a variety of life history 
traits. 
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